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We show that a clean multiband superconductor may display one or several phase transitions with 
increasing temperature from or to frustrated configurations of the relative phases of the supercon- 
ducting order parameters. These transitions may occur when more than two bands are involved in 
the formation of the superconducting phase and when the number of repulsive interband interac- 
tions is odd. These transitions are signalled by slope changes in the temperature dependence of the 
superconducting gaps. 
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The possibility of multiple bands contributing to the 
formation of a superconducting phase has been consid- 
ered in the case of transition metals [TH3] , supercon- 
ducting copper oxides [4] and magnesium diboride [5HZ]. 
More recently, sign-reversed two-band superconductivity 
has been proposed for the iron-based layered pnictides 
(HUH]. In the two-band case, the relative phase of the su- 
perconducting gap functions associated to each band is 
determined by the sign of the interband interaction, be- 
ing zero (tt) for attractive (repulsive) interband coupling. 
However, if more than two electronic bands have to be 
considered in the study of the superconducting phase, the 
relative phases are not uniquely defined by the signs of 
the interband interactions and in particular, frustration 
may occur if the number of repulsive interband interac- 
tions is odd. 

There is a close analogy between such frustrated multi- 
band superconductors and the well studied problem of 
frustrated Josephson junction arrays since the interband 
pairing may be regarded as an interband Josephson tun- 
nelling |10j . For example, it is known that a squared 
Josephson junction array with a 7r magnetic flux per pla- 
qucttc is frustrated with a degenerate ground state. The 
effect of the tt magnetic flux is the change of the sign 
of the Josephson coupling from positive to negative, i. 
e., a 7r-junction |llj . Such 7r-junctions are also present 
without magnetic flux as a consequence of d-wave pairing 
symmetry. In the case of a frustrated (odd number) loop 
of 7r-junctions, a spontaneous current will be present. 

Studies of frustrated Josephson junction arrays assume 
usually symmetric junctions. The case of multiband su- 
perconductivity is more closely analogous to the case of 
an array of asymmetric Josephson junctions [12] . When 
temperature is increased, strong modifications of the ra- 
tios of the interband Josephson tunnelling rates may oc- 
cur due to the relative changes of the superconducting 
gaps which are known to happen in multiband super- 
conductors [3 . We show in this manuscript that this 
may lead to one or two phase transitions with increasing 
temperature from or to frustrated configurations of the 
relative phases of the superconducting order parameters 
(which correspond to degenerate ground states which are 



chirally different) . In very recent works, the ground state 
of a three-band superconductor with repulsive interband 
interactions has been studied using a phenomenological 
Ginzburg -Landau approach |13j and a simplified BCS 
gap equation system where all intraband interactions are 
zero [13]. The effects reported in this manuscript are 
not observed neither in the temperature range where the 
Ginzburg-Landau approach is valid (since the gap func- 
tions have constant ratios in this temperature range) nor 
when all intraband interactions are zero. Furthermore, 
while the BCS gap equation does indicate the maxima, 
minima and saddle points of the free energy, it does not 
indicate which solution corresponds to the absolute min- 
imum of the free energy. 

In this paper, we address the case of a three-band su- 
perconductor with one or more repulsive interband cou- 
plings. This study can easily be generalized to any num- 
ber of superconducting bands. 

We adopt the Hamiltonian introduced by Suhl, 
Matthias, and Walker [3] to describe two-band super- 
conductors, generalized for n bands, 
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with = £fc — /i and where i, Vij, and /i are respectively 
the band index (i — 1, . . . ,n), the interaction constant 
and the chemical potential. The k sums in the inter- 
action term follow the usual BCS restrictions. S-wave 
symmetry is assumed throughout the paper for simplic- 
ity. The interaction constants are assumed to be positive 
(attractive interactions) for electrons in the same band 
(i — j), and one has Vij = Vji- 

In this paragraph, we follow the seminal discussion of 
Ref. [10] ■ We start by restricting the rt-band supercon- 
ducting system to the subspace constructed from the set 
of BCS states, 
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where A = (Ai, • • • , A„) and 3? = (</>i, • • • , </>„) are the 
absolute values of the superconducting gaps and the re- 
spective phases, and with the usual BCS definition of 
superconducting order parameters 
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with different values of A^ are orthogonal in the thermo- 
dynamic limit since their inner product gives a Dirac 5- 
function [10]. Defining the operators = J^k c ! 



which are diagonal in the previous basis in the thermo- 
dynamic limit, 
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one may write a simple expression for the eigenvalues of 
the Hamiltonian which is also diagonal in this basis, 
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where fid^il ) is the kinetic energy contribution of the 
respective band term in the Hamiltonian and has an elab- 
orate dependence on the superconducting parameters Aj. 
The other terms result from the intraband and interband 
pairing terms in the Hamiltonian. 

The phases associated with superconducting parame- 



ters, 



only affect the interband contribu- 



tions for the ground state energy. The minimization of 
this energy with respect to the phases <pi gives 
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These equations motivate a simple analogy with a 
system of n two-component classical spins governed 
by the Heisenberg Hamiltonian H = X)i>j ' Sj- 

The two-component classical spin is written as Si — 
Si(cos(j>i,sin<f>i). The energy of this system for a given 
set of magnitudes {Si} and angles {4>i} is Ej({4>i}) — 
<Aj) = J2i>j Jij cos((f>i - <pj) where 



J2i>i JijSiSj cos( 



Jij = JijSiSj. The last expression of this energy can 
also be interpreted as the energy of a triangular circuit 
of Josephson junctions with couplings Jy. The classi- 
cal ground state configuration is obtained from the min- 
imization of this energy with respect to the angles fa 
which leads to the condition £^ Jij sin(fa — fa:) = 
for each i, analogous to Eq. [6] with the correspondence 
Jij — > — 2V^|^>,||V>j|. Note that if a set of angles {fa} 
is a solution of the previous set of equations, then the 
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FIG. 1. (c) and (d): Plots of the Josephson energy Ej as 
a function of the phases fa and fa for respectively: (a) a 
frustrated array of three classical spins ; (b) a non-frustrated 
array. The Josephson couplings are indicated in (a) and (b). 
Two additional solutions of the minimization conditions (red 
spheres in the contour plot) are present in the frustrated case. 



set {— cf>i} also is. This implies that non-frustrated so- 



lutions must have 



= 0, ±7r. Other values cor- 



respond to frustrated configurations (degenerate ground 
states which are chirally different). 

We will now restrict our study to three bands (n = 3) 
but the arguments are easily generalized to any n. We 
impose for a matter of convenience that fa = (with 
no loss of generality). Solving Eqs. [6] we get several so- 
lutions corresponding to extreme or saddle points of the 
interband energy contribution, the non-frustrated solu- 
tions being ( fa, fa, fa) = (0,0,0), (0,7r, 0), (0,0, 7r), and 
(0, 7r, 7r) and the frustrated solutions being 



fa, fa) =± 0,cos 1 {a ),-sgn(^)cos 1 (a + ) 



where 
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and a = J12/J23 and b = JzijJiz- These frustrated 
solutions exist only if \a \ < 1. 

In Fig. [lja) and Fig. [ljb) we show a frustrated array of 
three classical spins and a non-frustrated array, respec- 
tively. In Fig. [TJc) and Fig. [ljd) , the respective plots of 
Ej as a function of the phases fa and fa as weu as the 
respective contour plots are displayed. The location of 
the solutions of the minimization equations is indicated 
by the red spheres in the contour plots. As J31 is var- 
ied from 1 to .2, the position of the frustrated solutions 
in Fig[TJc) continuously converges to the saddle point 
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is the renormalization of the interband coupling. 

The same reasoning can be followed at finite temper- 
ature. One must consider the free fermion entropy con- 
tribution for the free energy as well as the non-zero oc- 
cupation of quasiparticle states. One has 
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FIG. 2. The J31/J23 versus J12/J23 phase diagram of three 
classical spins when one of the interactions Jij is nega- 
tive. The non-frustrated regions have (<f>i, fa, fa) = (0,0,0), 
(0, 7r, 0), (0, 0, 7r), and (0, 7T, 7r). The boundaries of the regions 
of frustration correspond to second-order phase transitions. 



[fa, fa) = (0, 7r) [and to the equivalent location (0, — 7r)] 
and disappears when reaching this point [with (0, 7r) be- 
coming a local minimum] . This is a typical second-order 
phase transition. 

In Fig. [5J the J31/J23 versus J12/J23 phase diagram of 
the three classical spins system is displayed in the case 
when one of the interband Jij is negative. Note that 
in this case, the ratios of the couplings determine the 
signs of the couplings fully. The phase diagram of the 
three classical spins system when one of the interband 



and the interband pairing term has the same expression 
as that of Eq. [5j Since the other terms in the free en- 
ergy do not depend on the superconducting phases, the 
minimization with respect to the superconducting phases 
generates the same set of equations as for zero temper- 
ature (Eq. 6]). Minimizing the free energy with respect 
to the absolute values of the superconducting parameters 
[TOl H5| IT6] , one obtains a system of coupled gap equa- 
tions 
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with j = l,...,n, which, following the usual steps fl5| 
116) . can be written as 
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where E = y^£ 2 + A 2 , loq is the usual frequency cutoff, 
Jij is positive can be easily obtained from the former Nj(g) is the density of states of the j band and t = k B T. 



using the transformation 03 — > 03 + n which leads to 
^23 ~^23 and J31 — > — J31. For example, the sym- 
metrically frustrated case J12 = J23 = J31 = 1 has 
the frustrated solutions (fa, fa, fa) = ±(0, 27r/3, — 27r/3). 
These solutions using the previous transformation cor- 
respond to the frustrated case J 12 = 1 and J23 = 
J31 = —1 and (fa, fa, fa) = =p(0, 7r/3, 2ir/3). Recall- 
ing the relation with the multiband superconductor, one 
has J31/J23 -> (V31/V23) ■ (|*i|/|* 2 |) and J12/J23 -> 
(V12/V23) • (1*11/1*31), so that the location of the multi- 
band superconducting system in this phase diagram is 
determined not only by the interband pairings but also 
by the superconducting gaps. 

We have discussed in the previous paragraphs the min- 
imization with respect to the superconducting phases. 
Given a certain phase configuration, the remaining min- 
imization with respect to the absolute values of the su- 
perconducting order parameters gives the following result 
for n bands, 

A,fe = ^2 V kl> cos{fa - fa)u jk ,v jk , , (9) 
so that the effect of the superconducting phase differences 



We assume equal constant density of states for all bands, 
Nj(£) = Nj(0) = N, as a simplification. The differences 
in the density of states could also be absorbed in the 
couplings definition. 

The finite temperature minimization process with re- 
spect to the superconducting phases is the same as that 
described for zero temperature (with the same analogy to 
the classical spin system), but now the couplings Jij are 
temperature dependent. Let us consider the case where 
the interband couplings, Vy, are finite, but much smaller 
than the intraband pairings, Vu- The bands are indexed 
according to their uncoupled critical temperatures (crit- 
ical temperature associated with each band when all in- 
terband couplings are zero), starting at the highest T c 
and ending at the lowest. 

In the uncoupled situation, the superconducting gaps 
display sharp decreases at the respective critical temper- 
atures with temperature regions where only A3 is zero 
(between T c3 and T c2 ) and both A 3 and A 2 are zero 
(between T C 2 and T c ±). The modifications relatively to 
the uncoupled case, when the interband couplings Vij 
are small, are larger in these temperature ranges since 
the coupling of the gap equations implies that if Ai is 
non-zero at a local minima, then A2 and A3 are also 
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FIG. 3. (a) Superconducting gap functions and (b) super- 
conducting phases as function of the normalized temperature 
for a three-band superconductor with one repulsive and two 
attractive interband couplings. The vertical dotted lines de- 
limit the temperature range where frustration occurs. Slopes 
changes are clearly observed in the lower superconducting 
gap at the boundaries of this interval indicating second-order 
phase transitions, (c) Path followed (dashed orange curve) in 
the J31/J23 versus J12/J23 phase diagram by the correspond- 
ing three classical spin system. Parameters: V12 = —0.0045, 
V22 = .95, V23 = 0.016, V31 = 0.016 and V33 = 0.85 in units 
of V11. 



non-zero. So in the non-frustrated situation one expects 
for the three-band superconducting gaps an analogous 
temperature dependence to that typical of a two-band 
superconductor (with an extra gap function) [3]. 

However, the coupling depends strongly on the 
values of the superconducting parameters 14^1 and 
which, for weak interband interactions, follow approx- 
imately the behavior of the uncoupled gap functions, 
i.e., undergo sharp decreases around the respective un- 
coupled critical temperatures. Therefore one should ex- 
pect a fast increase in the absolute value of J31/J23 — > 
(V31/V23) • (l^il/l^l) around the uncoupled critical 
temperature T c2 and a similar increase of J12/J23 — > 
(^12/^3) • (l^i l/l^al) around the uncoupled critical tem- 
perature T C 3. Such variations may lead to crossings 
from a frustrated region of the phase diagram to a non- 
frustrated region or vice- versa. This is shown in Fig. |3|c) 
for a three-band superconductor which is initially non- 
frustrated and completely crosses the frustration region 
with increasing temperature. This system has weak 
(one repulsive and two attractive) interband couplings: 
V12 = -0.0045, V22 = -95, V 23 = 0.016, F31 = 0.016 



and V33 = 0.85 in units of V\\. The plot of the su- 
perconducting gap functions as function of the normal- 
ized temperature for this system is shown in Fig. |3ja). 
The associated superconducting phases are displayed in 
Fig. [3^b) (only one of the frustrated solutions is shown, 
the other being symmetric). Both graphs were obtained 
solving numerically the minimization conditions, Eqs. [6] 
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In Fig. [3]ja), slopes changes are observed in 
ower superconducting gap which reflect a second- 
order phase transition from or to a frustrated phase con- 
figuration as shown in Fig. |3jb). The vertical dotted 
lines indicate these transitions. Note that tiny slope 
changes also occur in the other superconducting gap 
curves. These changes are very small because the in- 
terband contribution to the respective gap function val- 
ues is much smaller than the intraband contribution. 
In Fig. |3jb), the superconducting phases do not reach 
the point {<j>\, 4>2, <fe) = T(0, 7r/3, 2tt/3) corresponding to 
maximum frustration. This is because the \Jij\ are never 
simultaneously equal to one as one increases temperature. 

The transitions to or from frustrated configurations 
are more difficult to occur if the interband pair tun- 
nellings is of the order of or stronger than the intra- 
band pairing, since the ratios of the superconducting 
gaps in this case have little temperature dependence 
leading to short paths in the phase diagram of Fig. [2] 
This can be illustrated considering the particular case 
when Vij — a 1+ ^~ 2 Vu, with |a| < 1, which has a non- 
frustrated solution A3 = a/S.2 — oj 2 Ai, which implies 
that the ratios J31/J23 and J12/J23 are constant and 
therefore the path in the phase diagram of Fig. [2] be- 
comes a single point. One should emphasize that even 
with weak interband couplings, in the case described 
above of one repulsive and two attractive interband cou- 
plings, if the three-band superconductor is initially in the 
(0, 7T, 0) non-frustrated region, no transitions will be ob- 
served with increasing temperature. 

As far as we know, there are at the present no experi- 
mental reports of the temperature-induced phase transi- 
tions reported in this manuscript. One should recall that 
necessary conditions for the existence of such effects are: 
i) more than two bands participate in the superconduct- 
ing state; ii) attractive intraband interactions; iii) weak 
interband pair tunnellings; iv) odd number of repulsive 
interband interactions. Furthermore, the second-order 
phase transitions to or from frustrated configurations are 
only clearly observed in temperature ranges where the in- 
terband contribution is dominant at least for one of the 
gap functions. Among the several examples of supercon- 
ductors where the possibility of multiple bands contribut- 
ing to the formation of a superconducting phase has been 
considered, the iron-based layered pnictides seem to be 
the more suitable candidates for the observation of these 
effects due to the complex band structure with more than 
two Fermi surfaces (hole and electron like). 

In conclusion, we have studied a three-band supercon- 
ductor considering frustration effects due to the existence 
of repulsive interband interactions. With increasing tern- 
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perature and in the case of small interband couplings, frustrated configurations of the superconducting phases 
one may have second-order phase transitions to or from which lead to slope changes in the temperature depen- 
dence of the superconducting gaps. 
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